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ABSTRACT 

We present correlations between 9 CO transitions (7 = 4-3 to 12-11) and beam-matched far-infrared (far- 
IR) luminosities (Lfir.*) among 167 local galaxies, using Herschel Spectral and Photometric Imaging Re- 
ceiveriCaFourier Transform Spectrometer (SPIRE; FTS) spectroscopic data and Photoconductor Array Cam¬ 
era and Spectrometer (PACS) photometry data. We adopt entire-galaxy FIR luminosities (Lfir.c) from the 
IRAS Revised Bright Galaxy Sample and correct to Lfir,* using PACS images to match the varying FTS beam 
sizes. All 9 correlations between L'^q and Lfir.a are essentially linear and tight (cr = 0.2-0.3c/ex dispersion), 
even for the highest transition, 7 = 12-11. This supports the notion that the star formation rate (SFR) is 
linearly correlated with the dense molecular gas (hhi 10"*“® cm“^). We divide the entire sample into three 
subsamples and hnd that smaller sample sizes can induce large differences in the correlation slopes. We also 
derive an average CO spectral line energy distribution (SLED) for the entire sample and discuss the implied 
average molecular gas properties for these local galaxies. We further extend our sample to high-z galaxies 
with literature CO(7 = 5-4) data from the literature as an example, including submillimeter galaxies (SMGs) 
and "normal" star-forming BzKs. BzKs have similar FIR/CO(5^) ratios as that of local galaxies, and follow 
well the locally-determined correlation, whereas SMG ratios fall around or slightly above the local correla¬ 
tion with large uncertainties. Finally, by including Galactic CO(7 = 10-9) data as well as very limited high-z 
CO(7 = 10-9) data, we verify that the CO(7 = 10-9)-FIR correlation successfully extends to Galactic young 
stellar objects, suggesting that linear correlations are valid over 15 orders of magnitude. 

Subject headings: galaxies: ISM — galaxies: starburst — infrared: ISM — ISM: molecules 


1. INTRODUCTION 

Molecular gas is the raw material from which stars are 
formed, hence its presence and physical conditions correlate 
with the star formation rate (SFR). A well-known correla¬ 
tion between gas surface density and SFR surface density is 
the Kennicutt-Schmidt law (K-S law): T,sfr = with 

slope N = 1.4 ±0.15 (Kennicutt 1998), where I^sfr is based 
on Hq, luminosity or far-infrared (far-IR) luminosity (Lfir, 
40-400/im) and is the sum of the atomic gas (from 
HI observations) and molecular gas (inferred from CO(l-O) 
line luminosity) surface density. This quantitative correlation 
serves as a fundamental input to most cosmological simula¬ 
tions (e.g. Springel & Hernquist 2003; Mac Low & Klessen 
2004; Krumholz & McKee 2005; ICrumholz & Thompson 
2007; Narayanan et al. 2008). However, it brushes over 
and/or oversimplihes the relationship between the individual 
gas components and star formation (SF). First, observations 
have shown that the presence and state of the atomic gas com¬ 
ponent has little impact on SF (e.g. Bigiel et al. 2008; Schruba 
2013; Liu & Gao 2012; Liu, Gao & Greve 2015). Second 
the correlation between the most widely used molecular gas 
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tracer, CO(l-O), and Lfir is found to be more complex than 
just a single power law: Gao & Solomon (2004) showed that 
the CO(1-0)-FIR slope changes from 1.27 to 1.73 with differ¬ 
ent sample selections, e.g. including only normal star-forming 
galaxies (SFGs: Ljr < 10** L©), or adding (ultra-)luminous IR 
galaxies (LIRGs: Ljr = lO'^-'^L©, ULIRGs: Ljr > lO'^L©), 
respectively. Greve et al. (2014, hereafter G14) showed simi¬ 
lar slope variations from 0.9 to 1.4 by analyzing a large num¬ 
ber of data from the literature. It is thus questionable that 
a single power-law relation exists between CO(l-O) and FIR 
and whether it can be applied to all galaxies even if it exists. 

CO(l-O) has a low critical density (nH2,m/ ~ 2 x 10^ cm“^), 
and thus traces the total amount of molecular gas. It pro¬ 
vides no physical insight into the higher density gas (hhj > 
10"* cm“^) that is known to be found at the formation sites of 
individual stars (Kennicutt & Evans 2012). The sensitivity of 
the CO (7,,—>7„-1) transition to denser/warmer gas increases 
with an increasing upper level (7„): for example, CO(10-9) 
has a critical density of nyi^ cnt ~ 10® cm ^ (Carilli & Walter 
2013, hereafter CW13). In contrast, the ground-7 of high 
dipole-moment molecules such as HCN and HCO^ already 
have nFi 2 ,crit 10®cm“^, thus probing the dense gas in cold 
phase. Studies based on HCN (Gao & Solomon 2004; Wu et 
al. 2005, 2010) and high-7 CS (Zhang et al. 2014) have found 
a unihed linear correlation between dense gas and the SFR, 
describing the scenario in which all dense gas above a certain 
density threshold has a similar efficiency/timescale to collapse 
and form new stars, thus is linearly determining SFRs. 

Until recently, published mid-to-high-7 (7„ > 4) CO data in 
local galaxies were still scarce due to the challenges of ob¬ 
serving at frequencies > 450GHz from the ground. Bayet et 
al. (2009) performed large velocity gradient (LVG) modeling 
with ground-based J,, < 1 CO to extrapolate 7u > 7 transitions 
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for several local galaxies. They found that CO versus total IR 
luminosity (Ljir, 8-lOOO^m) correlations have decreasing 
slopes for increasing Ju, similar to what modelings and simu¬ 
lations predicted (Narayanan et al. 2008; Juneau et al. 2009). 

More recently, G14 presented the Herschel Spectral and 
Photometric Imaging Receiver (SPIRE; Griffin et al. 2010) 
Fourier Transform Spectrometer (FTS, Naylor et al. 2010) 
high-7 CO data from the HerCULES program (Van der Werf 
et al. 2010; Rosenberg et al. 2015). They find that up to CO(7- 
6) the CO versus TIR or FIR relations are all roughly linear, 
while for higher-7u the correlations become sub-linear, similar 
to the results of Bayet et al. (2009). 

In this work, we use a large SPIRE FTS data set to sta¬ 
tistically determine the correlations between 9 L'qq (7 = 4- 
3 to 12-11) and beam-matched Lfir.*, and to determine which 
transitions are the best tracers of star formation. We also 
verify the slope variations against 7„ and test the validity of 
dense gas versus SF relation at high-z. The sample, data, and 
method are described in Section 2. The results and a discus¬ 
sion are given in Section 3. We adopt Hq = 73kms~^, Ua = 
0.73, Qm = 0.21. 

2. SAMPLE AND DATA 

The sample was selected from all public FTS observations 
in the Herschel Science Archive (HSA) of local galaxies that 
are in the IRAS Revised Bright Galaxy Sample (Sanders et 
al. 2003). A further selection criterion requiring the avail¬ 
ability of 70- 160/im band imaging data taken with Photo¬ 
conductor Array Camera and Spectrometer (PACS; Poglitsch 
et al. 2010) for galaxies partially resolved by the FTS beam 
(see Section 2.1) was imposed to enable beam-matching tech¬ 
niques to be used in the analysis. The final sample contains 
167 local galaxies (z < 0.064, di < 286Mpc), including 124 
(U)LIRGs and 43 SFGs (see Rosenberg et al. (2015) and Lu 
et al. (2014) for details of the HerCULES and GOALS pro¬ 
grams, respectively; details of the full data set are given in 
Liu et al. in preparation). 

2.1. FTS CO Measurements 

FTS has two bolometer arrays: the SLW bolometer ar¬ 
ray (446-989GHz) and the SSW bolometer array: (959- 
1543 GHz). Two observing modes are used: single-pointing 
for point-like sources (mostly with angular size less than FTS 
beam sizes, and di > 3QMpc), and mapping for nearby ex¬ 
tended sources (di < 3QMpc). In single-pointing observa¬ 
tions the central bolometer of each array is coaligned with the 
target, while off-axis bolometers point to the off-center sky. 
Mapping observations perform jiggling to scan extended tar¬ 
gets, and each individual bolometer will produce one spec¬ 
trum at its own pointing (R.A., deck). All mapping data 
and ~ 30% single-pointing data were reduced using SPIRE 
vl2 calibration products and Herschel Interactive Process¬ 
ing Environment (HIPE vl2.L0, Ott 2010) pipelines, with 
the remainder reduced using SPIRE vlO calibration prod- 
uctsH-pipelines. We note that measured line fluxes for the two 
sets of the ~ 30% single-pointing data show little difference 
(cr ^ 8%). For single-pointing data, we use central bolometers 
to extract CO lines. For mapping data, which are usually for 
spatially-extended galaxies, we select individual bolometers 
that are >half-beam-separated (Nyquist sampling) to extract 
CO lines. 

Note that FTS covers a wide range of frequencies, and its 


beam size varies accordingly ^ 43" at CO(4-3) to ^ 17" 
at CO(13-12), and is not a simple function of frequency. The 
beam size variation is well-calibrated for each central bolome¬ 
ter, but unavailable for off-axis bolometers. Thus we assign 
an additional 15% uncertainty for off-axis bolometers. 

To derive CO line flux, we use two different line profile fit¬ 
ting functions in HIPE: a Sine and a Sinc-convolved-Gaussian 
(SCG). A SCG, where a Gaussian line profile is convolved 
with the Sine instrument response, is appropriate for obser¬ 
vations of sources with broad/resolved lines (e.g. Zhao et al. 
2013): in such cases a Sine fit would underestimate the line 
flux by ^ 40%. In the more common case where the CO line 
is not obviously resolved, the SCG-derived fluxes are system¬ 
atically larger by < 20%. Given that a large fraction of single¬ 
pointing data have broad/resolved lines while mapping data 
do not, here we use SCG and Sine for single-pointing and 
mapping data, respectively, although we note that using only 
SCG or Sine for all data has no significant change to our cor¬ 
relation results. The CO flux uncertainties are derived from 
the rms in the baseline-subtracted spectra in the vicinity of 
each CO line rather than from formal line fitting uncertainties. 
Then the fluxes are converted to L'qq according to Solomon et 
al. (1992). 

2.2. Beam-matched FIR 

As a starting point we took the entire-galaxy FIR luminos¬ 
ity (LpiR.e) from Sanders et al. (2003). Given the mismatch 
between the dramatically varying FTS CO beam sizes and 
the physical extent of nearby galaxies (e.g. Galametz et al. 
2013) and some merging/interacting (U)LIRGs (e.g. Gao & 
Solomon 1999), we used a beam-scaling method based on 
the PACS imaging photometry ^ to scale entire-galaxy Lfir.c 
down to the local region that matches the CO beam size area 
(Tfir.^, for each CO line of each bolometer). PACS bands 
span the peak of the spectral energy distribution (SED) of 
typical FIR-luminous local galaxies, and thus provide a good 
proxy for FIR and are insensitive to dust temperatures (Tdust) 
By performing photometry with an aperture of CO beam 
size (LpAcs.i) an aperture of entire-galaxy (TpAcs.c); re¬ 
spectively, we can determine the FIR luminosity that is appro¬ 
priate to each CO beam size: Lfir,;. = FpACS,b/FpACS,e>^LFiR,e- 

In this way we calculate a value of Lfir.* for each CO line. 
This is essential for local galaxies that are partially resolved 
for FTS CO beam sizes (e.g. G14; Bayet et al. 2009; Yang et 
al. 2013; Zhang et al. 2014). 

3. RESULTS AND DISCUSSION 
3.1. Correlations between Lqq and Beam-matched Lfir 

Fig.l shows the 9 CO-FIR correlations of local galaxies 
and spatially resolved regions of nearby galaxies. FTS can ob¬ 
serve CO(4-3) to CO(13-12) except for galaxies at z > 0.032 
(mostly (U)LIRGs) whose CO(4-3) shift out of the FTS SLW 
waveband. Moreover, detections of CO(13-12) are sparse, 
and thus are not analyzed here. 

To derive the slope and intercept, we fit all data points (ex¬ 
cluding upper limits) with two linear fitting codes: an IDL 
least-squares fitting code based on MPFIT (Markwardt 2009), 
and an IDL Bayesian regression code LINMIX_ERR (Kelly 

* Eigure 5.18 of http://herschel.esac.esa.int/Docs/SPIRE/spire_handbook.pdf 

® All PACS data are updated to calibration 48 and post-processed with 
Scanamorphos vl9.0 (Roussel 2013). 

PACS bands: 70, 100, and 160/rm. Apeak ~ 290 p,mjTdust, where rjust ~ 
20—SOS' for typical local SFGs. 


High-J CO IR Correlations 
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Figure 1. Nine coirelations between L^q and ipiR.* among local galaxies (blue squares and green triangles) and spatially resolved regions of nearby galaxies 
(orange circles). The blue squares are from the HerCULES sample, and are the brightest (U)LIRGs. The green triangles are the other single-pointing data which 
are less bright, consisting of (U)LIRGs, nearby Seyfert nuclei and spiral nuclei. The orange circles are the mapping data, of which many are nearhy galaxy 
nuclear regions, but still have some off-nuclear regions at the highest transitions. The numbers of detections for each (sub)sample are shown at top left, next to the 
symbol legends. The error bar of each symbol legend is the mean uncertainties of the (sub)sample. The arrows are Scr upper limits. The best-fitting parameters: 
slope N, intercept A, Pearson correlation coefficient R and scatter a (dex) are shown at bottom right. The solid lines are the best-fit lines with free slopes, while 
dashed lines are the hest-fit lines with slopes fixed to one, which represent the average Tfir/Z-^q for each transition. In the last panel CO(12-l 1) we show the 
best-fit line of G14 as a blue dotted line for comparison. 


2007). The two methods give consistent results, thus we only 
show the latter results in Fig. 1 . 

CO(4-3) has the largest number of detections, but the least 
dispersions are seen in the CO(6-5) and CO(7-6) panels (both 
with tj = 0.\9dex). This conhrms the conclusion of Lu et al. 
(2014) that CO Ju ^ 6 — 1 transitions are the best SFR tracers 
for (U)LlRGs. 

For the entire sample, we find no evidence of signihcant de¬ 
crease in slopes with increasing N ^1.0-1.1 at ^ 4-6, 
and 0.96-1.0 at Ju 10- 12. Besides, the httings are con¬ 
sistent with the average distribution of upper limits for non¬ 
detections in each transition. 

Bayet et al. (2009) analyzed the IR beam corrections based 
on 850^m images and found correction factors similar to ours 
for galaxies we had in common. They found that IR beam 


corrections have only minor effects on the slopes (< 5%) and 
thus did not apply them, while in contrast we found that slopes 
without IR beam corrections are smaller by 10-15%. This 
partially explains the discrepancy between our results and 
their decreasing slopes, i.e. ~ 1 at ~ 3, then N ~ 0.8 
at Ju ~ 7, and hnally N = 0.53 ± 0.07 at Ju = 12. 

In addition, small number statistics is another key factor 
behind the discrepancies between different studies. G14 pre¬ 
sented the high-7 CO-FIR correlation with HerCULES FTS 
data (Van der Werf et al. 2010) and data the from literature 
of high-z submillimeter galaxies (SMGs). Their HerCULES 
sample contains 23 (U)LIRGs (excluding extended/merging 
galaxies). They evaluate IR beam correction with 870/im or 
350/fm maps, but the correction factors are small for their lo¬ 
cal sample (mostly (U)LIRGs). They found correlation slopes 
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1.0 for 1 < 7„ < 5, which then decrease to A^ = 0.87 ±0.05 
at Ju = 1 and then rapidly become N = 0.51 ± 0.11 at /„ = 12. 

For comparison, we divide our entire sample into three 
subsamples: a HerCULES subsample corresponding to 
the G14 local sample (mostly (U)LIRGs, excluding ex¬ 
tended/merging ones), all other single-pointing data (mixed 
SFGsh-(U)LIRGs), and all mapping data (nearby resolved 
SFGs, dominating the faint-end). Table. 1 lists the best- 
fit slope N, the intercept A, and the mean FIR/CO (A = 
log(LFiR/Fco)) for the entire sample and three subsamples. 
For the HerCULES subsample, similar to G14 and Bayet et 
al. (2009), we find decreasing slopes, e.g. N = 0.65 ± 0.31 at 
CO(12-ll). We overlay the significantly sub-linear CO(12- 
11)-FIR correlation of G14 {N = 0.51 ± 0.11) in the CO(12- 
11) panel of Fig. 1 for reference. 

Thus, with the largest high-/ CO data set available to date, 
we conclude that within the luminosity range shown in Fig.l, 
L'co ^ 9 closely and linearly follow LpiR and by exten¬ 

sion the SFRs, while all other 4 < 7,, < 12 CO-FIR correla¬ 
tions are tight and not far from linear. These results are not in 
conflict with previous studies. 

3.2. Average CO/FIR Spectral Line Energy Distribution 
(SLED) 

In Fig.2, we show a global FIR-normalized CO SLED, 
constructed from the products of the inversion of the best- 
fit FIR/CO normalization parameter for each transition (A, 
which is equal to Lfir/L^q, given in Table 1) (where 
Ju^ is included to make the unit same as integrated line flux 
Jy km ). Error bars are the dispersions of CO-FIR correla¬ 
tions (the cr shown in Fig. I). 

This CO SLED represents the average CO excitation con¬ 
ditions in local galaxies, and reveals at least two components: 
a low-excitation component peaking around J,, ^3-4, and a 
higher-excitation component peaking around 7„ ^ 8. 

We used RADEX (van der Tak et al. 2007) to construct 
two-component LVG models and performed least-x^ fitting to 
search for the best-fits of gas kinetic temperature Tpin and den¬ 
sity hhjj deriving 7i;„ ~ 90;!^“ K, log hhj ~ 3.0;!;o3 cm~^ for 
the lower-excitation component, and Tki„ > 200 K, log nH 2 
4.1;|;q 2 cm~^ for the higher-excitation component. 

In Fig.2, we also overlay the FIR-normalized CO SLEDs of 
a selection of galaxies that have at least six CO detections and 
span a wide Lfir range. (U)LIRGs/starbursts have stronger 
high-7 CO excitations and their normalized SLEDs are flat¬ 
ter at 7„ > 6, whereas normal SFGs with weaker excitations 
show a bump at 7„ ^ 4 in their normalized SLEDs. Surpris¬ 
ingly, the largest difference between normalized SLEDs of 
(U)LIRGs and SFGs is seen not at a high-7 (i.e. 7„ ~ 10) but 
at a low-7 (i.e. 7a ^ 4), where the CO/FIR ratios vary within 
^1.5 dex. Since the low-excitation component dominates the 
low-7 part of SLED, the large variation at a low-7 indicates 
that low-excitation gas is less correlated with FIR, while the 
high-excitation component is intrinsically better correlated 
with FIR (see also Lu et al. 2014). 

3.3. The CO(5^) - FIR Correlation and Extension to High-z 

To extend the high-7 CO-FIR correlation toward higher-z, 
we use the extensive compilation of high-z CO data in CW13 
that contains 60 CO 7„ > 4 detections with literature Lir". 

** The LpiR for nine galaxies are derived from IR SEDs, a further 14 from 
FIR-radio correlation (CW13; Bothwell et al. 2013), and the remaining 37 


1 - CO SLEDs Normalized by FIR 


Colored by log,gL^,p 



'iHn (low excitotion) (high excitolion) 

1^1 I I I 1^1 I I 


1 2 3 4 5 6 7 8 9 10 11 12 13 

J 

upper 


Figure 2. Average CO/FIR SLED (i.e. CO SLED normalized by FIR) for the 
entire sample, overlaid with several FIR-normalized CO SLEDs of individual 
galaxies spanning a wide range of Lfir. The black open squares are the values 
of the average CO/FIR (i.e. L^q/Lfir x 7^ ). The black solid curve is the best- 
fit two-component LVG model. The other dashed curves are the best-fit one- 
or two-component LVG models for individual galaxies. The colors indicate 
their Ig(LFiR). The two embedded panels are the least-x^ fitting results of 
low-excitation gas (left) and high-excitation gas (right) respectively. 



Figure 3. CO(5—4)-FIR correlation including high-z galaxies. The cyan 

diamonds are high-z galaxies from CW 13. The crossed red circles are z ~ 1.5 
BzKs from Daddi et al. (2015). The other data points and the best fit of the 
local sample are the same as the CO(5-4) panel of Fig. 1. 

They are the most extreme starbursts across all cosmic time, 
likely experiencing a short-term burst phase of SF resembling 
local ULIRGs, rather than the long-lasting mode in normal 
SFGs or the so-called ’’main-sequence” (MS) galaxies at high- 
z(e.g. CW13). 

coy„ > 4 observations toward high-z MS galaxies are still 
very rare: Daddi et al. (2015) presented the first CO(5^) 

scaled from single-band 850/im or 1200/im (CW13; G14). Where required 
we converted Lxir to Lfir t>y a factor of 1/1.3. 
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Table 1 

Best-fit parameters of Beam-matched CO-FIR correlations 


(Sub)sample 

CO(4-3) 

CO(5-4) 

CO(6-5) 

CO(7-6) 

CO(8-7) 

CO(9-8) 

CO(10-9) 

CO(ll-lO) 

CO(12-ll) 

Entire 

A =1.06 ±0.03 

1.07 ±0.03 

1.10±0.03 

1.03 ±0.04 

1.02 ±0.03 

1.01 ±0.03 

0.96 ±0.04 

1.00 ±0.05 

0.99 ±0.09 


A =1.49 ±0.24 

1.71 ±0.22 

1.79 ±0.24 

2.62 ±0.26 

2.82 ±0.27 

3.10±0.22 

3.67 ±0.25 

3.51 ±0.33 

3.83 ±0.67 


A =1.96 ±0.07 

2.27 ±0.07 

2.56 ±0.08 

2.86 ±0.07 

3.04 ±0.08 

3.20 ±0.09 

3.38±0.10 

3.56±0.11 

3.77±0.15 

HerCULES 

A =0.75 ±0.62 

1.08 ±0.43 

0.98 ±0.29 

0.84 ±0.35 

0.81 ±0.24 

0.81 ±0.26 

0.73 ±0.25 

0.76 ±0.26 

0.65 ±0.31 


A =4.46 ±0.57 

1.67±3.89 

2.78 ±2.61 

4.18±3.06 

4.68 ±2.04 

4.80±0.21 

5.56 ±0.02 

5.49 ±0.02 

6.58 ±0.00 


A =2.20 ±0.23 

2.40 ±0.21 

2.61 ±0.21 

2.85 ±0.21 

3.08 ±0.21 

3.28 ±0.22 

3.42 ±0.22 

3.65 ±0.23 

3.90 ±0.28 

Pointing 

A=1.00±0.10 

1.10±0.09 

1.06 ±0.09 

1.05 ±0.08 

1.09±0.11 

1.02±0.11 

0.93 ±0.09 

1.05±0.12 

1.12±0.18 

A=1.96±0.88 

1.43 ±0.79 

2.09 ±0.79 

2.51 ±0.65 

2.34 ±0.92 

3.02 ±0.84 

3.88 ±0.70 

3.12±0.91 

2.73 ± 1.36 


A=1.93±0.13 

2.30±0.12 

2.62±0.11 

2.89±0.11 

3.02±0.13 

3.18±0.16 

3.34±0.16 

3.51±0.18 

3.65 ±0.24 

Mapping 

A=1.17±0.10 

1.13±0.09 

1.11±0.10 

0.96 ±0.08 

0.94 ±0.08 

1.17±0.12 

1.12±0.18 

1.03±0.17 

0.66 ±0.42 

A =0.64 ±0.71 

1.29 ±0.67 

1.78 ±0.74 

3.08 ±0.56 

3.44 ±0.56 

2.14±0.77 

2.68±1.14 

3.40 ± 1.04 

5.99 ±0.00 


A =1.93 ±0.08 

2.22±0.10 

2.47 ±0.13 

2.82±0.12 

3.04±0.12 

3.19±0.12 

3.40±0.17 

3.58±0.18 

3.81 ±0.29 


Note. — The best-fit CO — FIR correlations for the 9 transitions. We fit the entire sample and three subsamples separately: G14 FlerCULES data, 
single-pointing data, and mapping data. Both x and y eiTors are considered. Upper limits are not considered in the fit. N is the best-fit slope and A is the 
best-fit intercept from free-slope fitting. A is the mean FIR/CO ratio (A = log(LFiR/L^Q)). 


detections in a sample of 4 z ~ 1.5 BzK-color-selected MS 
galaxies (BzKs). Thus combining all high-z > 4 data, 
CO(5^) has the second-most detections but covers the most 
diverse galaxy types. 

In Fig. 3, we show the CO(5^)-FIR correlation, including 
high-z galaxies as an example to illustrate the correlation be¬ 
tween dense gas and SFR, as CO(5-4) traces a hundred times 
denser molecular gas than CO(l-O). We show the linear fits 
from Fig. 1 only, as the inclusion of high-z SMGs and BzKs 
does not significantly change the results of the fit. The mean 
log(LFiR/LpQ( 5 _ 4 ^) = 2.27±O.O7L0 (Kkms~^ pc^)~^ is consis¬ 
tent with Daddi et al. (2015) considering a conversion factor 
of 1/1.3 from TIR to FIR. BzK galaxies fall on the local linear 
best-fit correlation (dashed line), whereas high-z SMGs are 
offset above by ^ 0.28 t/ex. Interestingly, using HCN(l-O) as 
a dense gas tracer, Gao et al. (2007) found higher FIR/HCN 
ratios in high-z galaxies (mostly QSOs/AGNs with only 5 
HCN detections) than those of local ULIRGs/SFGs. 

Given the large uncertainties in Lfir for a large fraction of 
SMGs compared to BzKs (from SED fitting), the amount 
of excess in FIR/CO(5^) (i.e. dense gas SF efficiency) in 
SMGs should still be treated with caution. And this excess is 
not significant enough to break down a linear form of dense 
gas versus SF relation, as these starbursts are the most extreme 
systems at all redshifts and as such do not represent the domi¬ 
nant mode of SF observed in MS galaxies (e.g. BzKs). Thus, 
better FIR measurements, e.g. via SEDs, and additional high- 
J CO observations in normal MS galaxies are needed before 
making solid conclusions at high-z. 

3.4. The CO(10-9) - FIR Correlation Extending to Galactic 
Young Stellar Objects (YSOs) 

In Eig.4 we show the correlation between an even 
denser/warmer gas tracer CO(10-9) and Lfir, which has the 
largest number of observations of sources spanning high-z 
galaxies (from CW13 and ALMA Partnership et al. (2015)) to 
galactic YSOs/protostars (from San Jose-Garcia et al. (2013)). 

ALMA Partnership et al. (2015) presented spatially re¬ 
solved ALMA CO(10-9) in a strongly lensed SMG SDP.81 
at z = 3.042, with demagnified Ltir = 5.1 x lO'^L©, compa¬ 
rable to the BzKs. We show the global SDP.81 as well as its 

For example, 850/rm-derived Lfir (using CW13 conversion factor) are 
~ 77-200% of radio-derived Lfir in Bothwell et al. (2013). Tensing correc¬ 
tions also introduce uncertainties. 



L'co [X km s"' pc’] 

Figure 4. CO(10-9)-FIR correlation, including galactic YSOs/protostars 

(San Jose-Garcia et al. 2013, red triangles), high-z SMGs/QSOs (CW13, cyan 
diamonds), and strongly lensed SMG SDR81 (global and two components; 
ALMA Partnership et al. 2015, red crossed circles). The 3 data points of 
SDP.81 are: the west component, the east component, and global, from left 
to right, respectively. The solid line is a free-slope fit for all sources. The 
dashed line is a fixed-slope A = 1 fit, indicating a mean log (Lfir/L^q^iq^^j) = 
3.30 ±0.09. 

east and west components in Pig.4. 

The best fit for all sources is A = 1.04 ± 0.01 and A = 
3.04±0.06 (see also San Jose-Garcia et al. 2013, who ob¬ 
tained the same N = 1.04 with only six galaxies); consid¬ 
ering YSOs/protostars alone we get N = 1.26 ±0.09, which 
is likely biased by the limited number of bright sources 
{Lhoi ~ 10^ Lq). a fixed-slope A = 1 fit indicates that a mean 
log(LFiR/L 0 Q(jg_gj) = 3.30±0.09 LQ(Kkms~^pc^)~^ is valid 
for all sources within a ^ 0.36dex. 

3.5. Conclusions 

We use Herschel SPIRE LTS observations of 167 local 
galaxies, including mapping data, to determine L!qq 7 = 4-3 
to 12 - 11 , and derive corresponding beam-matched Ljur i, for 
each CO line using PACS photometry. We find that these 
CO-LIR correlations are all linear, and that the non-linear re¬ 
sult reported in G14 can be attributed to the comparatively 
small number of galaxies in their sample. The overall lin¬ 
earity suggests that a universal CO/LIR SLED exists among 
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these galaxies, which further reveals two excitation states of 
gas. The high-excitation CO SLED peaks at 7 ^ 7, where the 
transitions are the best tracers of SFR (e.g. Lu et al. 2014). 

At high-z, however, SMGs/QSOs have elevated FIR/CO(5- 
4) ratios relative to those seen in local galaxies and contempo¬ 
rary BzKs: this is because they are the most extreme types of 
galaxies and hence are not representative of the typical SFG 
population. 

A tight, linear correlation between CO(10-9)-FIR is shown 
to hold for the majority of local galaxies, resolved sub-kpc re¬ 
gions, Galactic YSOs, and also high-z galaxies. These results 
strongly support a fundamental linear relationship between 
dense gas and the SFR, and also provide the local benchmark 
for probing the gas and SF in more types of galaxies (e.g. MS 
galaxies) at high-z. 


This work is supported by NSFC #11173059, #11390373, 
#11420101002, and CAS #XDB09000000. D. L. gratefully 
thanks T. Greve, Z. Zhang, R Papadopoulos, S. Madden, and 
R. Wu for constructive discussions, and K. Okumura and B. 
Altieri for helpful instructions on PACS. 
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